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We calculate the 0{ol 2 s (3q) corrections to the decay rate b — > ccs. For rea- 
sonable values of m c /mb this term is of the same order as both the one-loop 
^ ■ and 0{p? s log 2 my/ /mi) corrections to the decay rate. For m c /m]y = 0.3 the 

h ; 

O(«g/?o) corrections enhance the rate by ~ 18% . We also discuss the 0(a?Po) 
corrections to R T , the B semileptonic branching fraction and the charm mul- 
tiplicity. 
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I. INTRODUCTION 



The doubly-charmed decay mode of the B meson, B — > X c5s , has been the object of recent 
interest, since this mode makes a significant contribution to the inclusive B semileptonic 
branching fraction [ttlQ. Recently, the one-loop corrections to b — > ccs were calculated ||||] 
(see also ||) and found to be substantial, giving a ~ 22% enhancement to the tree level rate 
(for rh c = 0.30). This is significantly larger than the corresponding ~ 5% 0(a s ) correction 
to b — > cud decay. When combined with the additional radiative corrections, this brings 
the theoretical prediction for the semileptonic branching fraction into agreement with the 
experimental observation, within the theoretical uncertainties ||. 

Since the typical energy released in the decay, A = mj, — 2m c (neglecting the s quark 
mass), is much less than m&, one might expect the relevant scale for the perturbative correc- 
tions to b — > ccs to be significantly less than m&. Indeed, as stressed in Ref. f7|, the energy 
release in this process is so small that the assumption of local duality may not hold; it has 
been argued in Ref. || that deviations from duality would not show up at any finite order in 
the operator product expansion. However, even if the assumption of local duality does hold 
in this instance, this low scale would result in an even greater enhancement of this mode over 
the tree-level result. This is a higher order effect which requires a full two-loop calculation 
to address, which we have not attempted. However, in the approach of Brodsky, Lepage 
and Mackenzie (BLM) |§ useful information may be obtained by simply calculating the rif 
dependent piece of the order a 2 contribution to the decay. This determines the contribution 
of O(a 2 f3 ), where /3 = 11 — |n/. Since /3 is large, this term dominates the two loop result 
for many processes. The BLM scale /j^blm for the one-loop correction is defined as the scale 
at which the O(a 2 f3o) correction is absorbed in the one-loop correction. This approach has 
recently been used to estimate the two-loop corrections to semileptonic top, bottom and 
charm decays [pL0|, 1 1 1 . 



In this paper we calculate the O{a 2 (3o) correction to the decay b — > ccs. We will find that 
this term enhances the decay rate by almost as much as the one loop term, and is of the same 
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size as the 0(a 2 s log 2 mw/mb) correction. However, as we will discuss, the (9(a 2 /3o) term is 
not necessarily expected to dominate the remaining uncomputed two-loop corrections. 

In Section 2 we compute the O{a 2 s f3o) corrections to the mode b — > ctV t . This contribu- 
tion arises from strong interaction corrections to the be vertex, and the result can be related 
to a piece of the b — *> ccs correction by taking m T = m c . These corrections are interesting 
in their own right as they give O(a 2 /3o) corrections to the ratio 

_ r(6 -> X c rv T ) 
Rr = T(b - X c eu e ) • (L1) 

We will find that the two-loop corrections to this ratio are under control. In Section 3 we 

calculate the 0(ct 2 /?o) corrections to the cs vertex. We give our conclusions in Section 4. 

II. 0(a| fa) CORRECTIONS TO b -► ctV t 

The rates for B — > X c ev e and B — > X c tv t may be written as power series in a s and 
h-Q,CT>/ m b [0- The leading order result in l/m b reproduces the parton model, while to 
0(l/ml) two unknown nonperturbative parameters, A and Ai, arise. The ratio R T defined 
in Eq. Ql.l] ) provides a potential constraint on these parameters, although the uncertainty 



in the measurement is currently too large for these constraints to be useful [15-18]. As in 



the case with massless leptons JTTJ] the O(a 2 (3 ) corrections to this process are quite large; 



however, these corrections largely cancel in the ratio R T . 

We write the semitauonic decay of a 6-quark in terms of the quark pole masses m;, and 



m c as 

r(6 - crv r ) = ^frW (m c , m T ) ( 1 + ^IrW (m c , m T ) {2.1) 

where rh c = m c /mb, rh T = m T /mb, fio = 11 — ^n/ is the QCD /3-function, and rif is the 
number of light quark flavors running through the vacuum polarization loops. The ellipsis 
denote terms O(a^) and higher. The one-loop correction r^ 1 - ) (m c ,m T ) is given in Ref. ||T4| . 



To compute the O(a 2 fto) term we follow the work of Smith and Voloshin [10] and 



compute the O (a s ) rate with a finite gluon mass, T^ 1 \m g ,m c , rh T ). The 0(a 2 Po) correction 
in the true theory, T^ 2 (m c ,m T ), can be found from this rate by performing the weighted 
integral 

rf(m c1 rh T ) = -Po^^ r ^ fr«K,m c ,m T ) - -^^r«(m c , J (2.2) 
p Arc Jo m 2 g \ m l g + m l h J 

where a^\mb) is the strong coupling defined in the ^/-scheme of Ref. B, and is related to 



the coupling a s defined in the MS scheme by 

a ( y\^) = a s ^) + ^^-P + ... . (2.3) 

We have obtained a lengthy analytic expression for dT(m g )/dq 2 , where y/q 1 is the invariant 
mass of the lepton pair, which we have integrated numerically over q 2 and m g to obtain ri 2 \ 
Since the results are very sensitive to rh c and rh T , we have chosen to follow the approach of 



Refs. [18] and express these ratios as a power series in l/m B : 

Ac = ! « £ _Af 1 _5£)_ll( 1 _^) + + (2 .4) 

m B m B V m B J m B \ m B J 2m B m D \ m B J 

A A 2 A 

= 0.372 - 0.628 0.628^- + 1.16— 

m B m B m B 

m T ( A A 2 Ai \ 
m T = — 1 + + 



Tn B y m B m 2 B 2m 2 B J 

A A 2 A 

= 0.334 + 0.334 + 0.334^- - 0.167— i- , 

m B m B m B 

where we have defined the spin-averaged meson masses 

4 2m D V ; 

m B + 3m B * , Ai , Q1QH 

m B = = m b + A h . . . ~ 5313 MeV. 

4 2m B 

To the order in which we are working we can just use the leading term in our perturbative 
calculation. We find 



T(B - X c tV t ) = \V bc \ 2 [0.082] 



192tt 
'a s (m b y 2 

7T 



A / a s (mh) 

1 - 1.94 1.29 " K ' 



m B \ 7r 

1.28 (^) (3 + O(l/m%,a s /m B ,a 



(2.6) 



For completeness, we also give the result for m c = 0.3 and m b = 4.80 GeV, 



T(B -> X c rV 7 



. 2 (j^mfc 
1 bc| 192tt 3 



[0.114] 



1 - 1.39 



' a s {m b ) 



+ 0(l/ 



71 



2 2 

m bl a s 



1.58 ( ^A) p 



7T 



(2.7) 



As is the case for b — > cev e decays, the O(a 2 p ) corrections in Eqs. (|2.6|) and fl2.7|) are 
quite large, corresponding to a low BLM scale for this process. However, these corrections 



largely drop out of the ratio R T . Combining Eq. Q2.6Q with the results of [pTlJI , we find 



R T = 0.224 
~ 0.224 



1 + 0.24^ + 0.15 (^A 



71 



71 



\ 2 T 

Po- 0.29 + 0(l/ 

J m B v 



m 2 B ,a 2 s 



(2.8) 



A 



1 + 0.017 + 0.007 - 0.29— + O { 1/m 2 , a 2 

m B v 



where we have taken a s (m b ) = 0.23 in the second line. The perturbation series appears well 
behaved, and the corresponding BLM scale for R T is /zblm = 0.29 m&. 



III. O{a 2 s p ) CORRECTIONS TO b -» ccs 

Neglecting the s quark mass[], we write the width for b — ► ccs decays (where the final 
state includes an arbitrary number of gluons and light quarks) as 

64vr 3 



r(6 - cc.) = (m c ) ( 1 + ^ r w (™ c) + ' r( 2 ) (m c ) + . . . ] (3.1) 



where 

T (0) (x) = Vl - 4x 2 (l - 14x 2 - 2x 4 - 12x 6 ) + 24x 4 (l - x 4 ) In ( 1 + ^ ~ Ax * \ (3 2 ) 
v y v y \1 — v 1 — 4x 2 / 

is the tree-level result; r^°^(0.30) = 0.196. The complete one-loop corrections may be ob- 
tained from Refs. [|,(§; for m c = 0.30 one obtains rW(0.3) = 2.99. Taking a s (m b ) = 0.23, 
this corresponds to a 22% enhancement of the rate over the tree level result. 



1 Since m s ~ Aqcd; we will treat terms of order m 2 to be of the same size as terms of order Aq CD , 
which we are neglecting. 
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Since the four-quark operators responsible for nonleptonic b decays run in the effective 
theory below mw the 0(a 2 ) contributions to the decay are more complicated than for 
semileptonic decays. We write the 0(a 2 ) contribution to the decay rate as 

r (2) (m c ) = ci In 2 + c 2 (m c ) In + c 3 (m c )/? + c 4 (m c ) (3.3) 
m b m b 

where c% — 4 jl9| . The subleading log contribution C2(m c ) was calculated in Ref. J|; for 



m c = 0.30 these authors find C2 = 3.34. 

Clearly the requirement that C3 dominates the two-loop correction, implicit in the BLM 
approach, will not hold in this process, since the non-vacuum polarization terms c\ and C2 
are enhanced by powers of Xnmw/ j m b . Separating these terms out, we may instead hope 
that C3 dominates over C4 due to the factor of (3q. However, even this assumption may 
not hold. Voloshin || has shown that in the limit in which the charm is produced nearly 
at rest C4 receives a large enhancement. For b — > cud, in this limit the analog of C4 is of 
order tt 2 , whereas for b — > ccs the Coloumb exchange graphs between the two slowly-moving 
charmed quarks give a contribution to C4 of order 7r 4 . While it is not known whether this 
enhancement is relevant for the physical value of rh c , it indicates that the O(a 2 /3o) terms 
need not dominate over the 0(a 2 ) terms. Nevertheless, as a first step towards understanding 
the size of the two-loop corrections to this process, we may calculate C3. 

While the complete series of leading and subleading logs has been summed to all or- 



ders ||19|j4|, we cannot consistently use these results since we are not summing all terms of 
O(ag\og n ~ 2 (m b /m w )P ). However, as was stressed in Ref. [0, \n(m w / m b ) ps 2.8 is not a 
large number, and the leading log expansion does not seem to work well for nonleptonic b 
decays. For example, for b — > cud decay the subleading O (a 2 \a(mw / n^b)) term is 2/3 the 
size of the leading O (a 2 \n 2 (mw/m b )^J term. Therefore, we choose to work consistently to 
O (a 2 ) and discard the rest of the leading and subleading log terms. The neglected terms of 
O(a^) and above are likely to be much smaller than the uncomputed 0(a 2 ) corrections. 

The calculation of C3 is simplified due to the fact that the graphs factorize into the 
contribution from the upper be vertex and the contribution from the lower cs vertex. The 



upper vertex contribution can be simply obtained from the corrections to b — > ctV t (by 
making the substitution m T — > m c ), while the contribution from the lower vertex require an 
additional calculation. 

For the lower vertex corrections, the kinematic structure of the phase space allows us 
to express the integrals over the momenta of the c quark, s quark, and gluon in terms of 
the spectral density of the charged V-A current (the imaginary part of the charged current 
vacuum polarization), 

rfoLK) = / M„Ml ImIV(g 2 ) dr 2 (P b ;p c , q) dq\ (3.4) 

where is the contribution from the be line and ImlLj^g 2 ) is the imaginary part of the 
vacuum polarization. The tensor structure of the vacuum polarization can be decomposed 
into a transverse and a longitudinal contribution, 

ImIV(g 2 ) = ^PM 2 ) + L, ~ PM 2 ). (3.5) 

Since the functions Pi{q 2 ) and Pt{q 2 ) depend only on the scalar q 2 , the integration over 
dT 2 (Pb] Pci q) can be carried out analytically with a simple computation, 



1 lower I m 9 ) ~ 7 
rn b 



K m b -m c ) -q [m b + m C/ 



Pi(q 2 ) + 



(m 2 b + m 2 c ) 2 + q\m 2 h + m\ - 2q 2 )] P t (q 2 )} ^A(m 2 , mj, g 2 )-^ (3.6) 

where A(x, y, z) = x 2 +y 2 + z 2 — 2xy — 2xz — 2yz. The resulting expression is quite lengthy and 
we do not present it here. The functions Pt{q 2 ) and Pi{q 2 ) have been previously calculated 
(for a massless gluon) to 0(a s ) in the context of QCD sum rules pl| . 

It is then a simple matter to integrate numerically the resulting expression over q 2 and 
m 2 to obtain c 3 as a function of m c . At the "reference point" m c = 0.3, we find c 3 (0.3) = 3.7. 
Using a s (m&) = 0.23 and /3q = 9, this corresponds to an 18% correction to the tree level 
result, almost as large as the one loop correction. The values of c^{m c ) for a range of values 
of m c are given in Table | and plotted in Figs. [I] and ^] along with the separate contributions 
from the upper and lower vertices. 



m c 


1 v >(m c ) 


c?,{rh c ) 


— ^-^1 v ' [m c j 


/ a 3 (m b ) \ 2 n n ( ~ \ 


0.20 


0.99 


0.83 


0.07 


0.04 


U.zo 


i.y i 


O 1 1 

Z.Li 


n 1 a 
U. 14 


n in 

U. 1U 


0.30 


2.97 


3.67 


0.22 


0.18 


0.35 


4.25 


5.81 


0.31 


0.28 


0.40 


5.85 


8.89 


0.43 


0.43 



TABLE I. Numerical values of the one and partial two loop corrections r^and C3 for b — > ccs 
decay. In the last two columns we have taken a s (mt,) = 0.23 and (3q = 9. 
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FIG. 1. Contributions to (dashed lines) and C3 (solid lines) from (a) the renormalization 



of the be vertex, (b) the renormalization of the cs vertex, and (c) the sum of (a) and (b). 
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FIG. 2. T^ 1 ' (dashed line) and C3 (solid line) as functions of m c (expanded view of Figure 1 (c)). 

It is useful to compare these results with the leading and subleading log corrections to 
;p( 2 ). F or jyi = 0.3 these are 



41n 2 ^ + 3.341n ^) f^V ^ (3.5 + 1.0)/% 



(3.7) 



where we have removed a factor of /3o = 9 to allow comparison with the second column in 
Table U For rh c = 0.3 the a 2 /?o term is roughly the same size as the 0(a 2 s ) leading log 
correction, and a factor of three greater than the 0{a 2 s ) subleading log. 

Note that the 0(a s ) and O(a 2 s f3o) corrections to the cs vertex are positive for all values 
of m c , while the corrections to the be vertex are negative. The one loop corrections cancel 
at rh c w 0.14, while the O(a 2 /3o) corrections cancel at a slightly higher value of rh c . In 
this situation the BLM scale /xblm is not physically relevant: at the point where the one 
loop corrections to the vertices cancel, /^blm is singular, whereas at the point where the 
O(a 2 /3 ) contributions cancel /^blm = w&. In this region the BLM scale for the decay width 
is unrelated to the BLM scales that would be obtained for the upper and lower vertices 
individually, and does not reflect the average momentum of the gluons in the diagrams. 
Therefore we prefer simply to present our results as a contribution to the 0(a 2 ) correction 



to the decay rate. 

We also note that since the leading order phase space function r^ ^ (rh c ) is very sensitive 
to the b and c quark masses, there is a large uncertainty in the total b — ► tics width simply 
due to the uncertainty in the b and c quark masses. Since m b and m c are related via the 
1 / itlq expansion to the corresponding hadron masses, this sensitivity is really an additional 
hidden source of I/tuq corrections, just as in the semileptonic decay width. This is made 
clear if we adopt the approach of the previous section and write m c as a series in 1/mg. In 
this case, the large sensitivity to m c results in a I/tub correction which is as large as the 
leading order term, 

A 



G 2 m 5 

T(B -> X c , s ) = \V bc V cs \ 2 -^f [0.057] 



9.7— + O (I /mi, a. 
m B \ 



(3.8) 



Of course, one could argue that this result is misleading because we are expanding about 
the extreme value m c = 0.37. Nevertheless, the large l/m b correction shows the sensitivity 
of the width to the quark masses. Working instead with pole masses and keeping mj, fixed, 
varying m c between 0.27 and 0.32 results in a factor of two change in the total rate. 

It is straightforward to find the a 2 ^o term for the decay b — > cud from computations of 
the charmed semileptonic decay [|ll|] and from the results for R e + e - [^0|. For m c = 0.3 this 
gives 



O 2 m 5 

r(b^dad) = \v he v ud \ 2 [o.52] 



64tt 3 



1-0.67 



a s {m h ) ^ ^ 2 mw_ ( a s {m b ) 



71 



+7.17 In 



m w ( a s (m b ) 



2 



7T 



m b 



71 



i.ii (^) A + o(^) 



(3.9) 



Combining this with the results of the present work, we find the ratio of the partial widths 
for m c = 0.3, 

2 



r(6 



CCS) 



F(b — > cud) 



0.376 



\v cs 


2 


v ud 





1 + 3.66 sV b> + 4.80/5o [ V ' 



71 



71 



(3.10) 



-3.83 In 



m w ( a s (m b ) 



m b 



+ . 



7i 



The a^Po correction enhances the tree-level ratio by 22%. 
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Similarly, the O(a 2 /3o) enhancement of T(b — *> ccs) will decrease the semileptonic branch- 
ing fraction and increase the charm multiplicity (n c ). Combining the result for b — > ccs with 
the O(a 2 j3 ) corrections to the other modes, we find for m c = 0.3, an O(a 2 j3 ) correction 
shift to the semileptonic branching fraction of 



Since we are simply illustrating the effect of the O(a 2 /3o) terms on these observables, we do 
not include the remaining perturbative corrections or contributions from rare decay modes 
in these expressions. 



We have computed the O (ct 2 Po) contributions to the rate of the nonleptonic decay b — > 
ccs at the parton level. While these corrections do not dominate in any formal limit of 
the theory, they are a well-defined subset of the complete two-loop corrections. When the 
perturbation series is expressed in terms of a s (m&), the O(a 2 (3o) corrections are of the same 
order as both the one-loop corrections and the leading log corrections. For rh c = 0.3 they 
provide an additional reduction of ~ 1% in the semileptonic branching fraction, and increase 
the charm multiplicity (n c ) by ~ 0.04. 

These corrections are sufficiently large to cast doubt on the applicability of perturbative 
QCD to this decay mode. Since there is so little phase space, this is not unexpected. These 
corrections are in addition to the large 0(a 2 ) corrections suggested by Voloshin ||, as well 
as the large implicit l/m^ c corrections due to the uncertainties in the c and b masses. 




(3.11) 



The corresponding shift to the charm multiplicity (n c ) is 




(3.12) 



IV. CONCLUSIONS 
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